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being restricted to it, exemplified for the case of hot QED. The nonperturba- 
tive resummation obtained in a <I>-derivable approach is gauge independent, 
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I. INTRODUCTION 

One of the main issues of the ongoing heavy-ion program is the investigation of decon- 
fined hadronic matter. Despite asymptotic freedom of the strong interaction, this quark- 
gluon plasma is characterized by a large coupling in the regimes of physical interest, so 
nonperturbative approaches are required to describe this many-particle system reliably. 

Recently, within the framework of an equilibrium description of the QCD plasma, the 
calculation of thermodynamical quantities by resumming hard-thermal-loop (HTL) propa- 
gators was proposed. In |I|| and |Q] , nonperturbative expressions were given for the thermo- 
dynamical potential which, however, are detracted from inconsistencies: The expressions do 
not reproduce, in the weak-coupling limit, the perturbative results at leading order 0{g'^) 
and they are, with uncompensated medium-dependent divergences, not satisfactory from a 
formal point of view. These problems were claimed to be solved only in a fully resummed 
calculation. In Ref. 0, instead, a consistent approximation of the entropy was derived 
from HTL propagators. This approach in principle resolves the problem of a leading order 
thermodynamical resummation since, up to an integration constant, the thermodynamical 
potential can be reconstructed from the entropy ^. 

Nonetheless, a direct calculation of the thermodynamical potential, the quantity contain- 
ing the full thermodynamical information, is instructive, in particular since an approach to 
derive macroscopic properties, which are sensitive to hard momenta, from HTL propagators 
(as a soft-momentum approximation) is, a priori, far from obvious. On the other hand, 
being explicitly gauge independent and respecting for arbitrary momenta the fundamental 
sum rules resulting from the commutator relations, the HTL propagators are a favorable 
basis for calculating physical quantities from approximate dressed Green's functions, yield- 
ing (mostly) analytical nonperturbative results. These points are addressed in the present 
note where, starting from the Luttinger-Ward representation of the thermodynamical po- 
tential 0, the hot QED plasma is studied. This case is particularly simple to analyze but 
at the same time also representative for systems with an equal HTL structure as, e.g., the 



quark-gluon plasma at finite temperature and density. 

II. THE APPROXIMATION SCHEME 

As an exact relation, the thermodynamical potential can be expressed in terms of fully 
dressed Green's functions by the (generalized) Luttinger-Ward representation 



n = -Q[D]-n[s] + <^[D,s]. (1) 

fi is a functional of the photon propagator D or the electron propagator S, which are related 
by Dyson's equation to the respective self-energies. The boson part, e.g., is defined by 

i n[D] = ^ Tr [In(-D-i) + DU] , D-' = D,' - U , 

where the trace is taken over the four-momentum and the Lorentz structure, while in the 
analog fermion part — 17[S'], with S^^ = Sq'^—T,, the spinor indices are traced. The functional 
^[D, S] given by all two-particle irreducible bubble graphs with exact propagators ('dressed 
skeletons') is related to the self-energies by 

Consequently, the fundamental stationarity of the thermodynamical potential upon variation 
of the self-energies [§] is fulfilled by the representation (|l]). It is emphasized that, on account 
of its stationarity, Q is formally gauge independent, although the propagators are not. 

Using the projectors P^^ = -K^.KJK'^ and Pj, = g^,u - K^K^/K^ - P^^, where K = 
[K{Ku) — uK'^]/[{Kuy — K'^Y^'^ and u is the medium four-velocity, the inverse photon 
propagator is decomposed into the transverse (T) and the longitudinal (L) part as well as 
the covariant gauge-fixing term, 

i=T,L ^ 

with /S.Q^ = K^ = kQ — k^. Introducing the fermionic 'projectors' V±{K) = | (^ ± I^) on 
the particle and the hole excitations (the index denotes the ratio of chirality to helicity), the 
electron propagator can be written in a similar way as 



i=± 
In terms of the scalar propagators Aj, with the degeneracy factors (It = d — 1, di = 1 and 
d± = {d+ l)/2, the Q parts of eq. (|l|) read 



I ^[D] = l$l^^d^ [in(-Ari) + A,n,] - ln(-A 
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f!|Sl = X E <*. Iln(-Ar') + A.Ao ' " d^ l.i(-A„-') J , (3) 

where the subtractive contribution of the ghost fields, which otherwise decouple, is included 
in fi[D], and the integral-sums, continued to (i = 3 — 2e spatial dimensions, 
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.. ^ tr ' ^^'^ " ^ (2^)' 



run over either bosonic or fermionic Matsubara frequencies /cq- 

Commencing from this exact approach, selfconsistent ('symmetry conserving') approxi- 
mations can be derived : an approximation of the functional $ in the scheme (^^ yields 
approximate self-energies and an expression for the thermodynamical potential which is still 
stationary. In particular, the loop expansion of $ can be truncated at a certain order, which 
is, in terms of the perturbative expansion in free Green's functions, equivalent to a partial 
resummation avoiding the problem of double counting of diagrams. The leading-loop order 
of the ^-derivable approximation is given diagrammatically byQ 



which implies the relation 



2 2 ^ ^ 



among the selfconsistent solutions of the coupled Dyson equations and $ . 



^It is noted that the (vacuum) counter terms required to render the self-energies finite can be im- 
plemented in the ^-derivable scheme by appropriate counter loops in <I>, leaving fi itself unchanged. 



In the following, the selfconsistent one-loop self-energies are approximated by the HTL 
self-energies 11* and S*. Although the HTL self-energies remain a reasonable approximation 
even for hard momenta [|l^,|ri|], it is not evident a priori that a consistent approximation 
of Q can be formulated in terms of such approximated quantities which are derived for soft 
momenta much smaller than the temperature T, whereas thermodynamics is sensitive to 
the momentum scale T. In fact, the HTL approximation of the self-energies undermines the 
self consistency of the ^-derivable approach since the relation analogous to @ is not ful- 
filled. Nevertheless, as shown in the following, it is possible to consistently resum the HTL 
contributions to the thermodynamical potential, yielding a nonperturbative 'continuation' 
of the 0{e^) perturbative result. In Sec. HI, the HTL contributions to the Vt parts of the 
thermodynamical potential are calculated directly from the HTL self-energies since, e.g., 
replacing H'^ by H* is correct up to terms beyond the order 0{e^) under consideration. The 
corresponding $ contribution is then given in Sec. IV to complete the resummed approxima- 
tion of the thermodynamical potential which is, expressed by HTL self-energies, explicitly 
gauge invariant. 



III. THE n CONTRIBUTIONS 



The HTL self-energies of the photon and the electron are given by [|12 



H*=M2 + II, H*=-2II, Ii{k,,k) = Ml^ 
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M^ 
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can be considered as asymptotic masses (squared) of the transverse photon and the electron 
particle excitation, respectively, since their dispersion relations approach mass shells for 
momenta k ^ eT. The longitudinal photon (plasmon) mode and the hole (plasmino) 



excitation, on the other hand, possess a vanishing spectral strength when approaching the 
hght cone exponentially fast for k ^ eT. 

Besides the $ part, eq. (|^) evaluated with HTL propagators yields the contributions 
to the thermodynamical potential which are, up to gauge group factors, the same as for 
non-Abelian gauge theories with HTL self-energies with a structure like (^. In Refs. [|I],0|, 
however, where the QCD plasma was studied in the framework of the leading order HTL 
perturbation theory, only the ln(— Aj~^) terms of (|^) were considered to contribute to the 
thermodynamical potential. In terms of the bare-propagator expansion, this incomplete 
analysis amounts to a miscounting of graphs, so in the weak-coupling expansion the lead- 
ing order perturbative contribution to the thermodynamical potential is not reproduced 
correctly in either the case of hot QCD with vanishing chemical potential [0] or for the 
degenerate quark-gluon plasma at T = 0. For the QED plasma under consideration, the 
corresponding boson (6) and fermion (/) contributions, marked here by the index A, follow 
from the expressions given in |I|] by replacing the asymptotic gluon and quark masses by 

the expressions (^, schematically 

1 r. M'^ 1 

^h = 2 ?C {^^ In(-AV') + In(-Ar') - H-^o'} = -^- + finite terms , 

nij =^d± {in(-AV^) + ln(-A^_-^) - In(-A^i)} = finite terms . (7) 

The imperfect cancelation of the thermal divergences in Q\ ;, is a second, formal indication 
of missing contributions in the approaches [|l],§ . 

The complete analysis of the HTL-resummed thermodynamical potential has to keep 
track of the remaining terms, indexed B in the following, of eq. (|^). Using complex contour 
integration, the Matsubara sums can be calculated to yield a quasiparticle part stemming 
from the pole uji{k) of the propagators, and a Landau-damping contribution arising from 
the discontinuity of the HTL self-energies below the light cone. 






rk /7(M 

+ / — (1 + 2n{u;/T))^, , ^, = Disc(A^/-,) , 
JO Zir 

where n = in^j is either the Bose or the negative of the Fermi distribution function, 
and fi can be H^^ or A,^^ Subtracting and adding the appropriately infrared- regularized 
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asymptotic integrands analog to the technique apphed in [|l|, the integrals over the spatial 
momenta are split into a finite contribution and the dimensionally regularized subtraction 
term. Taking the limit e — > 0, the individual contributions are 
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with the angles \1/t,l = Disc(A^2,n^L)5 ^± = Disc(At|_Ao ^), and the subtraction terms 

4(A;2 + M2)\\ 
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It turns out in eq. (|^) that the thermal divergences of the fermion integral-sum cancel, as 
for ^*Af in (0). The boson part, on the other hand, contains a temperature dependent 
term ~ M^/e which is twice as large in magnitude as its counterpart in (|^). These thermal 
divergences of the A and B contribution have to cancel the corresponding terms in the 
nTL approximation of $, which is calculated in the following section, to yield a well-defined 
resummation of the thermodynamical potential. 



IV. THE * CONTRIBUTION 

In this section, the remaining $ contribution is calculated by evaluating the two-loop 
functional, given the self-energies and the Green's functions. | Within the leading loop 
approximation, $ is related to the self-energies by equation (H) and can thus be expressed 
in a general form by $'' = — itTrZ}''n" + | (1 — t) TrS'^S", independently of t. In the 
framework of the HTL approximation, however, the consistency relation (^ is violated since 
the traces are dominated by hard momenta. The resulting ambiguity is even of order O(e^), 
namely 
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Hence, the HTL contribution to $ cannot be obtained by naively replacing the leading 
loop quantities in $", as given above, by their HTL approximation. Instead, the HTL 
contribution $* can be conceived by analyzing how the C(e^) discrepancy arises. Denoting 
the photon momentum in Q by i^ and the fermion momenta by Qi^2, this diagram (with 
bare propagators for the C(e^) contribution) can be represented as a double integral-sum 
over an expression with a numerator A^ = K^ — Q\ ^ Q\- Closing the external legs of 
the boson self-energy in the HTL approximation amounts to neglecting the term K^ in N . 
Tracing over the negative of the fermion HTL self-energy, on the other hand, neglects one 
of the Q^ terms in A^. Thus, all terms are accounted for twice in the sum over all three 
possibilities to approximate one of the momenta as soft. Accordingly, the representation 

4 2 ^ ' 

reproduces $'' perturbatively to order 0{e^). As shown in the following, this representation 
indeed leads to a well-defined resummed approximation of the thermodynamical potential. 



^In the approach [^,^, instead, the functional properties of <1> in the leading loop approximation 
are used first to derive the entropy, which is then evaluated with the HTL propagators. It is noted, 
however, that in this framework the selfconsistency relation only holds approximately, see below. 



It is first emphasized that in the complete expression resulting from the HTL approxi- 
mation of eq. (H) and (pi), which can be written in a compact form as 



2 






Tr 



Zi 



f^ghost , (10) 



with the individual contributions given by (^, (||), all temperature dependent divergences 
cancel. Moreover, the perturbative limit of the thermodynamical potential is reproduced 
by the representation (plQl). Separating the free contributions, e.g., for the boson part by 
ln(— D~^) = ln(— Do^^)+ln(l— DoH), and using the expansion ln(l—a;)+x/2 = — s/2+C(a;^), 
the leading order correction to the interaction-free limit f^o = —{dT + | 2(i±) |q VT^ is 



1 1 T2 

- TrDon^ + - Tr^oS* = - 



^l = -- TtDoU* + - Tr^oS* = — (M^ + MJ) = fif^ . (11) 



As in the HTL calculations P,Q of the QCD entropy, the leading order term originates 
entirely from the behavior of the thermodynamically relevant excitations at the hard mo- 
mentum scale T, which a ;?ostenon justifies the present approach. However, in contrast to 
the entropy calculations |^,^ where the results are manifestly ultraviolet-finite, the cance- 
lation of the thermal divergences is in the present approach directly related to the fact that 
the perturbative result is reproduced. It is emphasized that this aspect makes the repre- 
sentation (p!0D of Q* unique; any linear combination, apart from (P), of closed self-energy 
diagrams for the ^-contribution would result in either uncompensated thermal divergences 
or an incorrect perturbative limit. 

The next-to-leading order term of the perturbative result Qp^^^, on the other hand, cannot 
be expected to be reproduced in the present approximation: In an equivalent approach for 
the scalar f^^^^ theory, where the complete leading loop Luttinger-Ward resummation of the 
thermodynamical potential can be derived |13[, the corresponding 0{g^) correction takes 
its correct value only after the resummation of the self-energy, while the expressions (^) are 
calculated with bare propagators. Accordingly, the contribution of order 0{e^), which arises 
from the static longitudinal parts of (0), 
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is found to underestimate the perturbative result fi^f^ = — (e^T^/3)^/^T/(127r) by a factor 
of 1/4. As observed in [|l],^ for QCD, the next-to-leading order term of ^^b agrees with 
the perturbative result but is overcompensated by the Q*^ ^ contribution, which is to be 
interpreted by a systematic next-to-leading order calculation. 

V. SUMMARY 

The generalized Luttinger-Ward representation of the thermodynamical potential is a 
suitable framework to derive consistently resummed approximations for relativistic gauge 
theories, with the propagators approximated by their HTL contributions. While the formal- 
ism is particularly simple to analyze in the case of the hot QED plasma, which is exemplified 
here, the application to other systems (including the QCD plasma) with the same HTL 
structure is evident. As a direct result of the HTL approximation, the resummed thermo- 
dynamical potential (|10|) is gauge independent. All medium-dependent divergences cancel, 
hence the approximation is explicitly renormalization-scale independent. The resummed 
expression (p!0| ) enjoys the anticipated behavior of a nonperturbative approximation. As 
shown in Figure 1, it yields a smooth extrapolation to the large-coupling regime, where it 
is inclosed by the perturbative results which are known to fluctuate with increasing order. 
In the weak-coupling limit, on the other hand, the perturbative result is recovered. This 
demonstrates that to leading order the resummed thermodynamical potential can entirely 
be expressed in terms of HTL propagators. 

The applications of the formalism to the hot and dense quark-gluon plasma are straight- 
forward and promising, in particular with regard to the suggestion [jl4| to extrapolate finite- 



temperature lattice data to non-zero chemical potential. The reliability of the leading order 
HTL resummation of the thermodynamical potential in the large coupling regime, however, 
remains to be justified by a systematic next-to- leading order calculation. 
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FIG. 1. The resummed pressure p* = —i}*/V (solid line), as a function of the coupling and 
scaled by the interaction- free limit po = (2 + |4) ^T^, compared to the perturbative result in 
leading order (dashed line) and in next-to-leading order (dash-dotted line). 
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